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Abstract 

We consider the dynamics of a harmonic crystal in d dimensions with n com¬ 
ponents, d,n > 1 . The initial date is a random function with finite mean density 
of the energy which also satisfies a Rosenblatt- or Ibragimov-Linnik-type mixing con¬ 
dition. The random function converges to different space-homogeneous processes as 
Xd ±00 , with the distributions p± . We study the distribution pt of the solution 
at time t G IR. The main result is the convergence of yt to a Gaussian translation- 
invariant measure as t —> 00 . The proof is based on the long time asymptotics of the 
Green function and on Bernstein’s ‘room-corridor’ argument. The application to the 
case of the Gibbs measures pL± = g± with two different temperatures T± is given. 
Limiting mean energy current density is — (0,..., 0, C(T+—T_)) with some positive 
constant (7 > 0 what corresponds to Second Law. 
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1 Introduction 


The paper concerns the problems related to the derivation of Fourier’s Law for harmonic 
crystals p|. We have started in |0-|^ the analysis of the long time convergence to the 
equilibrium distribution for partial differential equations of hyperbolic type in IR'’^ and for 
a harmonic crystal. Here we continue the analysis and prove Second Law for the crystal: 
the energy current is directed from high temperature to low temperature. Similar results 
have been established in Ba for d = 1, and we extend the results to all d > 1. The 
case d > 1 appears very different from d = 1 because of much more complicated properties 
of oscillatory integrals. This is why we combine here the methods from [O with new ideas. 


Namely, we develop our “cutoff’ strategy from |]^ which more carefully exploits the mixing 
condition in Fourier space. This approach allows us to cover all d > 1. 

We assume that the initial state Yq{x) of the crystal is a random element of the Hilbert 
space Ha of real sequences, see Dehnition 0 below. The distribution of 10 ( 3 ;) is a probabil¬ 
ity measure /xo of mean zero satisfying conditions S1-S3 below. In particular, the distribu¬ 
tion of Yo{x) converges to distinct translation-invariant measures /i± as Xd —>■ ±C )0 . Given 
t G IR, denote by the probability measure that gives the distribution of the solution 
Y{x, t) to dynamical equations with the random initial state Yq . We study the asymptotics 
of /it as f ±cx). 

Our main result gives the (weak) convergence of the measures /it on the Hilbert space 
Ha with a < —d/2 to a limit measure /ioo 


^ /^OO) t 


CX), 


( 1 . 1 ) 


which is a translation-invariant Gaussian measure on Ha ■ A similar convergence result 
holds for t ——CX) since our system is time-reversible. We construct generic examples of 
harmonic crystals and random initial data satisfying all assumptions imposed. The explicit 
formulas for the covariance of the measure /ioo are given in ( |2.14| )- (|2.23|) . We derive the 
expression for the limit mean energy current joo • The ergodicity and mixing of the limit 


measures /ioo follow by the same arguments as in |]TT 


We apply ( |1.1|) to the case of the Gibbs measures /i± = g± with two distinct temperatures 
T± > 0 (we adjust the dehnition of the Gibbs measures g± in Section 3). The measures g± 
satisfy all our assumptions, and the weak convergence gt goo follows from our results. We 
apply formula for the limit energy current joo = (ji,, ■ ■ ■, to the case of Gibbs measures 
g± and deduce that 

joo = -(0,...,0,G(T+-T_)), G>0. 

This corresponds to Second Law. 

For d = 1 similar problem have been analyzed also in [IT, IT]. The authors considered 
the hnite simple lattice of size L with the viscosity, in contact with two heat baths at 
temperatures T± . The convergence of the covariance is proved in the limit t > cx), and 
then L ^ 00 . The result is close to ours: the limit energy current is non zero and ~ AT 
which corresponds to the superconductivity p|. However, the space decay of the limit 


position-momentum covariance in [|lg is exponential which differs from the power decay 
in our problem and in [Q (see Remark |3.2| iii)). 

For case d > 1 the convergence 0 has been obtained for the hrst time in []^ for 
initial measures which are absolutely continuous with respect to the canonical Gaussian 
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measure. We cover more general class of initial measures with the mixing condition and 
do not assume the absolute continuity. For the hrst time the mixing condition has been 
introduced by R.Dobrushin and Yu.Suhov for the ideal gas 0. The condition substitutes 
(quasi-) ergodic hypothesis in the proof of the convergence to the equilibrium distribution, 
and plays the key role in our Bernstein-type approach. Developing this approach, we have 
proved the convergence for the wave and Klein-Gordon equations with translation-invariant 
initial measures ii p!4| . In IQ we have extended the results to the wave equation with the 
two-temperature initial measures. The present paper develops our previous results BH. 
where the harmonic crystal has been considered for all d in the case of translation invariant 
initial measure. Here we extend the results to the two-temperature initial measures. 


We outline our main result and strategy of proof. Consider a discrete subgroup T of IR'^, 
which is isomorphic to 2'^. We may assume T = after a suitable change of coordinates. 
A lattice in IR*^ is the set of the points of the form rx{x) = x-I-^a , where x G ^ IR”^, 

A = 1, ...,A. The points of the lattice represent the equilibrium positions of the atoms 
(molecules, ions,...) of the crystal. Denote by rx{x,t) the positions of the atoms in the 
dynamics. Then the dynamics of the displacements rx{x,t) — rx{x) is governed by the 
equations of type 


u{x,t) = - V{x -y)u{y,t), x G Z*^, 
j/eZ'* 

u\t=o = Uo{x), u\t=o = Vo{x). 


( 1 . 2 ) 


Here u{x, t) = (mi(x, t),..., Un{x, t)),uo = (mqi, • • •, won) Wo = (wi, • • •, ^'or^) G IR" , n = Ad; 
V{x) is the real interaction (or force) matrix, (Vki{x)'^ , k, I = l,...,n. Similar equations 
were considered in m, pTj pA|, p]] . Below we consider the system of type (|1.2|) with an 
arbitrary n = 1, 2 ,.... 

Denote Y{t) = {Y°{t),Y^{t)) = {u{-,t),u{-,t)), Yq = = («o(-)Wo(-)) • Then 


( p..2|) takes the form of an evolution equation 

Y{t) = AY{t), t G IR; Y(0) = Yq- 
Formally, this is the Hamiltonian system since 

AY = jA \\y = JVH(Y), ./=(_°iJ 


(1,3) 


(1,4) 


Here V is a convolution operator with the matrix kernel V and H is the Hamiltonian 
functional 

HiY) ■= ^{v, v) + ^(Vu, u), Y = {u, v), (1.5) 

where the kinetic energy is given by -{v,v) = - ^ |'?^(2^)P and the potential energy by 
-{Vu,u) = X X! {y{x — y)u{y),u{x)^ , stands for the real scalar product in the 


Euclidean space IR"^. 
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We assume that the initial correlation functions 

Ql\x,y):=E(Y^{x)®Y^{y)), x,ye 7 l.^, ( 1 . 6 ) 

have the form 

Qoix,y) = qo\x-y,Xd,yd), i,j = 0, 1 . ( 1 . 7 ) 

Here x = {xi,... ,Xd) = (x,Xd), y = {yi,..., yd) = {y, yd) e • Moreover, we assume that 

lim qo\z,yd + Zd,yd) = q±{z), z = {z,Zd) G (1.8) 

yd^±oo 

Here q±{z) are the correlation functions of some translation-invariant measures y± with 
zero mean value in Tia ■ The measure /tq is not translation-invariant if gj? 7^ . 

Next, we assume that the initial mean “energy” density is uniformly bounded: 

eo{x) := E[\uo{x)\‘^ + 1^0(2^)^] = tr( 5 o°(a:,x) -1- trQ“(x,a:) < Cq < 00, x E Z'^. ( 1 . 9 ) 

Finally, it is assumed that the measure /xq satishes a mixing condition of a Rosenblatt- or 
Ibragimov-Linnik type, which means that 


Yo{x) and Yo{y) are asymptotically independent as \x — y\ ^ 00 . 


( 1 . 10 ) 


To prove the convergence ( |1.1| ) we follow the strategy of There are three steps: 

I. The family of measures pt, f > 0 , is weakly compact in Ha , a < —d/2 . 

II. The correlation functions converge to a limit. 


Q7{x,y)= / [Y\x) ^Y\y)) pt{dY) ^ Q^^{x,y), t ^ 00 . 


( 1 . 11 ) 


III. The characteristic functionals converge to a Gaussian one, 

Pti^) := J ptidY) exp{-^Qoo(^, ^)}, t co. (1.12) 

Here T = (T°, T^) E V ■= D ® D ^ D := Go(Z‘^) (g) IR” , where Go(Z'^) denotes a space 

of real sequences with hnite support, (R, T) = I] S (W(x), T®(x)) , and Qoo is the 

*=o.WeZ'^ 

quadratic form with the matrix kernel (Q^J^(x, y))ij=o,i, 

E E (1.13) 

x,y&Z‘^ 


Below the brackets (■, ■) denote also the Hermitian scalar product in the Hilbert spaces 
L‘^{T'^) ® IR"" or its different extensions. 

For the proof of I - III we develop our cutting strategy from |H] combined with some 
techniques from [|l|. Each method and separately is not sufficient since here the 
measures pt are not translation-invariant and we consider all d > 1. To prove II we split 
Ql^{x,y) into even, odd components and the remainder as in [|l|. The even component 
corresponds to the translation-invariant initial measure and is analyzed by the method of 

0 


for all d > 1. On the other hand, the odd component is missing in |11| and it requires 
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a novel idea since its Fourier transform contains the Cauchy Principal Value which is more 
singular than measures corresponding to the even component. The singularity was studied 
in [m for the case d = 1. However, similar detailed analysis for d > 1 seems to be impossible 
due to the bifurcations of the critical points. 

Let us outline our method. We rewrite (|1.11|) in the equivalent form 


^ Qoo(d',T), t oo, 


(1.14) 


for 4/ G : by dehnition of T>^ , the Fourier transform T [6) vanishes in a neighborhood of 
a “critical set” C <ZT'^. The set C includes all points 6 & T'^ with a degenerate Hessian of 
, where ujI{9) are the eigenvalues of the matrix V{9) = . Also the set 

C includes the points 9 E either with uJk{9) = 0 , or Vg^uJk{9) = 0 or with non-smooth 
ujk{9). The cutting of the critical set C is possible by two key observations: i) mesC = 0 
and ii) the correlation quadratic form is continuous in P due to the mixing condition. The 
continuity follows from the space decay of correlation functions by well-known Shur’s lemma. 
The systematic application of the Shur lemma allows us to extend (|1.14|) from T G to 
all T G P by condition E6. 

Similarly, we hrst prove the property III for T G T>^ and then extend it to all T G P . For 
T G we use a variant of the S.N. Bernstein ‘room-corridor’ technique (cf. |]I[ for d = 1). 
We develop our variant of the S.N. Bernstein technique which we have introduced in |0- 


g], [Q in the context of the Klein-Gordon and wave equations and in for the harmonic 
crystal with d > 1 in the case of translation-invariant initial measures. For T G T>^ we have 
= E exp {i{Y{t),'^)). We rewrite, (V(t), T) = (V(0), $(-, t)), where <F(a;, t) andean 
be represented as an oscillatory integral. For <F(x, t) we get the uniform bounds (| 8 . 6 |) , (|8.7|) . 
These bounds follow by the stationary phase method because <F(x, 0) = T(a:) G and 
hence, T( 6 *) vanishes in all points 9 E C with degenerate Hessian of the phase function. 
The bounds roughly speaking imply the following representation: 




Y. Uy) 

y&Bt _ 

vW 


t ^ oo, 


(1.15) 


where Bt stands for the ball {y E Z'’* : \y\ < ct} and \Bt\ is its volume. Now (|1.12|) 
follows from ( |1.15|) by the Lindeberg Central Limit Theorem since To(|/i), To(|/ 2 ) are almost 
independent for large \yi — y2\ by mixing condition ( p..l 0 |) . 

Let us comment on our conditions concerning the interaction matrix V (a:). We assume 
the conditions E1-E4 below which in a similar form appear also in ^ |T^. El means the 
exponential space-decay of the interaction in the crystal. E2 resp. E3 means that the 
potential energy is real resp. nonnegative. E4 eliminates the discrete part of the spectrum 
and ensures that mes C = 0. We also introduce a new simple condition E5 for the case 
n > 1 which eliminates the discrete part of the spectrum for the covariance dynamics. It 
can be considerably weakened to the condition E5’ from Remark Hi). For example, the 
condition E5’ holds for the canonical Gaussian measures which are considered in [|T^. The 


conditions E4, E5 hold for almost all functions V{x) satisfying E1-E3 as shown in [pT 


Furthermore, we do not require that uJk{9) 7 ^ 0, as in [^j: note that a;(0) = 0 for the elastic 
lattice ( 2.26 ) in the case m = 0 . Our results hold whenever mes {6* G : ijJk{9) = 0} = 0 . 
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To cover this case we impose the condition E6 which is similar to the condition iii) from [13 


p.l71]. E6 holds for the elastic lattice (|2.26|) if either n > 3 or m > 0 . 

The main result of the paper is stated in Section 2 (see Theorem A). Section 3 concerns 
the application to Gibbs measures. In Section 4 we give bounds for the initial covariance. 
The compactness (Property I) is established in Section 5, convergence (p..ll|) in Sections 
6-7, and convergence (|1.12| ) in Section 8. In Section 9 we check the Lindeberg condition for 
convergence to a Gaussian limit. Appendix is concerned with a dynamics and covariance in 
Fourier space. 


2 Main results 

Let us describe our results more precisely. 


2.1 Dynamics 

We assume that the initial date Yq belongs to the phase space 1-La , « ^ 21^, dehned below. 


Definition 2.1 Ha is the Hilbert space of pairs Y = {u{x),v{x)) of -valued functions 
of X ^ endowed with the norm 

ll^lla^ E (|M(a:)|^ +|n(a;)p)(l +< cx). (2.1) 

We impose the following conditions E1-E6 on the matrix V. 

El There exist constants C,a > 0 such that \Vki{z)\ < /c,/Gn:={l,...,n}, z E 

7L^. 

Let us denote by V{Q) := (Vki{0)) where Vki{6) = yki{z)e^^^, 6 G T^, and T'^ 

\ / k,len -jjd 

z^£. 

denotes the d-torus IR'^/27rZ^. 

E2 V is real and symmetric, i.e. Vik{—z) = 14/(^) G IR, k,l En, z E'^^. 

The condition implies that V{6) is a real-analytic Hermitian matrix-function in 6 E T'^ . 

E3 The matrix V(9) is non-negative dehnite for each 9 E T'^. 

The condition means that the Eqn (|1.2|) is hyperbolic like wave and Klein-Gordon Eqns 
considered in Let us dehne the Hermitian non-negative dehnite matrix 

n{9) := {V{9)f^^ > 0 ( 2 . 2 ) 


with the eigenvalues /jJk{9) >0, k E n, which are called dispersion relations. For each 
9 E T'^ the Hermitian matrix fl{9) has the diagonal form in the basis of the orthogonal 
eigenvectors {ek{9) : k En} : 




n{9) = B{9) 

V 


0 

0 


0 

0 


\ 

B*{9), 


^ni.9) y 


(2.3) 
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where B{6) is a unitary matrix and B*{6) denotes its adjoint. It is well known that 
the functions ujk{0) and B{9) are real-analytic outside the set of the ‘crossing’ points 0*: 
0Jk{9*) = for some I ^ k. However, generally the functions are not smooth at the 

crossing points if 0Jk{9) ^ 0Ji{9). Therefore, we need the following lemma which is proved 
in B Appendix]. 

Lemma 2.2 Let the conditions El, E2 hold. Then there exists a closed subset C* C T'^ 
such that i) the Lebesgue measure of C* is zero: 

mesC* = 0. (2.4) 

a) For any point Q & \C^ there exists a neighborhood 0{Q) such that each dispersion 
relation uik{9) and the matrix B{9) can be chosen as the real-analytic functions in 0{Q) . 
Hi) It is possible to enumerate the eigenvalues uik{G) so that we have in the whole open set 

T'^ \ a 


(^) — • • • — (^)) a^ri+l (^) — • • • — ^r2 (^)) • • • ) (^) — • • • — (^)) (^•^) 

if a ^ u, I < r„,ru < Vs+i := n. (2.6) 

iv) Let us define Ila{9) for 9 E T'^ \ C^ as the orthogonal projection Ila{9) : IR"" — E„{9) 
onto the eigenspace E„{9) C IR" generated by the eigenvectors ek{9), k G {r„_i,r„] . Then 
no-(6'), a = 1, ..., s -1-1 , is a real-analytic function of 9 E T'^ \C^ . 


Below we suggest that u!k{9) denote the local real-analytic functions from Lemma | 2 . 2 | ii) 
Our next condition is the following: 

E4 Dki9)^0, Vfcen, where Dki9) := det 9 eT^\C,. 


Let us denote Cq := {9 E T'^ : detR( 6 ') = 0} and Ck := {9 E T‘^\C^, : Dk{9) = 0}, 
k = 1 ,... ,n . The following lemma has also been proved in [|TTi Appendix]. 


Lemma 2.3 Let the conditions E1-E4 hold. Then mes = 0 , A; = 0,1, ...,n. 


Our last conditions on V are the following: 

E5 For each k ^ I the identity ujk{9) — uii{9) = const_ , 9 E T^ does not hold with 
const _ 7 ^ 0 , and the identity ujk{9) -\- a;;( 6 *) = const+ does not hold with const + 7 ^ 0 . 

E6 ||R“^( 6 ')|| G L^{T'^) in the case when Co 7 ^ 0. 

This condition holds if Co = 0. 


The following Proposition is proved in [^, p.150], |^, p.l28] (see also Appendix). 


Proposition 2.4 Let El and E2 hold, and a G IR . Then 

i) for any Yq E Ha there exists a unigue solution Y{t) E CfWi^Ha) to the Cauchy problem 

B)- 

a) The operator U{t) : Yq i—> Y{t) is continuous in Ha ■ 


6 






2.2 Convergence to statistical equilibrium 

Let (n, S,P) be a probability space with expectation E and B{T-Ca) denote the Borel a- 
algebra in Ha ■ We assnme that Yq = Yq^uj, ■) in (|1.3|) is a measnrable random fnnction with 
valnes in {Ha, B{Ha)) ■ In other words, for each a: G the map uj i—>• Yq{ijj,x) is a mea¬ 
snrable map > IR^"' with respect to the (completed) a-algebras S and i3(lR^"'). Then 
Y{t) = U{t)Yo is again a measnrable random fnnction with valnes in {Ha,B{Ha)) owing 
to Proposition |^. We denote by /io(dLo) a Borel probability measure on Ha giving the 
distribntion of the Yq . Withont loss of generality, we assnme (ff, E,P) = {Ha, B{Ha), fJ^o) 
and Yq{u,x) = a;(a;) for /io(dci;) -almost all uj G Ha and each x G . 

Definition 2.5 /r* is a Borel probability measure in Ha which gives the distribution of 
Y{t) : 


pit{B) = /io(P(-t)P), VP G B{Ha), t G IR. 


(2.7) 


Onr main goal is to derive the convergence of the measnres /r* as f —*• cxo . We establish 
the weak convergence of fit in the Hilbert spaces Ha with a < —d/2 : 

'hi 

Rt —^ Roo as t ^ oo, (2.8) 


where /Too is a limit measnre on the space Ha , ol < —d/2. This means the convergence 

J f{Y)fn{dY) J f{Y)fiao{dY), t ^ cx), (2.9) 

for any bonnded continnons fnnctional / on Ha ■ 


Definition 2.6 The correlation functions of the measure Ht are defined by 

Q\\x,y) = E(Y\xH)®Y^{yH)), i,j = 0,1, x,y (2.10) 

if the expectations in the RHS are finite. Here Y^{x,t) are the components of the random 
solution Y{t) = (y'°(-, t), y'^(-, t)) . 


For a probability measure fi on Ha we denote by fi the characteristic functional (Fourier 
transform) 

/i(T) = J exp(i(F,T))/r(dF), T G P. 

A measure pi is called Gaussian (of zero mean) if its characteristic functional has the form 

/i(T) = exp{-iQ(T,T)}, T G P, 

where Q is a real nonnegative quadratic form in P. A measure pi is called translation- 
invariant if pi{ThB) = /i(P), P G B{Ha), G Z*^, where ThY{x) = Y{x — h) , x G Z'^. 
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2.3 Mixing condition 

Let 0(r) denote the set of all pairs of subsets .4., C at distance dist {A, B) >r and let 
(t{A) be a cr-algebra in generated by Y{x) with x & A. Dehne the Ibraginiov-Linnik 
mixing coefficient of a probability measure /xq on Ha by (cf |^, Dehnition 17.2.2]) 


V7(r) 


sup 

(A,B)£0{r) 


sup 

A e a{A),B e a{B) 
> 0 


lho(4 n B) — /xo(4)/io(j7)| 
ho(-B) 


( 2 . 11 ) 


Definition 2.7 The measure po satisfies strong, uniform Ibragimov-Linnik mixing condi¬ 
tion if 

ip{r) ^0 as r —>• cx). (2-12) 

Below, we specify the rate of decay of ip (see condition S3). 


2.4 Statistical conditions and results 

We assume that the initial measure po satishes the following conditions S0-S3: 

50 /io has zero expectation value, EYq{x) = 0 , x G 

51 po has correlation functions of the form (|1.7|) with condition (|1.8|) . 

52 po has a hnite mean energy density, i.e. Eqn (|1.9|) holds. 

53 po satishes the strong uniform Ibragimov-Linnik mixing condition with 

H-OO 


ip = J r'^ ^p^^‘^{r)dr < oo. 


(2.13) 


Introduce the correlation matrix of the limit measure poo ■ It is translation-invariant 
Q^{x,y) = = {qZi^ - 2/)),,^.= 0,1 ' 


In the Fourier transform we have locally outside the critical set C* (see Lemma |2.2|) 

&) = B{e)MZ{9)B*{e), 1,J = 0,1, (2.15) 

where B{6) is the smooth unitary matrix from Lemma p.2|, ii) and M^{6) is n x n-matrix 
with the smooth entries 


M^d{0)ki = Xki 
Here we set (see 


f)l ) 

(B-(«)(Ay)«(«)B(«))„ + isgn(^(9))(B'(9)(M„-)«(9)S(9)) 

D) 


kl 


. (2.16) 


Xkl — 


1 if k,l e {ra_i,ra], a = 1,..., s-b 1, 

0 otherwise 


(2.17) 









with ro := 0, r^+i := n , and 


<(») - 5(q+(«) + c’(«)q*(«)C'*(«)), 
M„-(e) := i(c(«)q-(«)-q“(«)(?•(«)). 


with q"*" := 


1 

2 


(?+ + ?-), q' 



g_) and 


C{e) := 


0 n-\e) 

-n{e) 0 


c*{e) 


0 -n{e) \ 

n-\e) 0 ) ’ 


(2.18) 

(2.19) 


( 2 . 20 ) 


where C* denotes a Hermitian conjugate matrix to the matrix C . The local representation 
( p.l5| ) can be expressed globally as the sum: 


fc(9) = +iioo(9). 


( 2 . 21 ) 


where 


s + l 


(&)”(«) := i:n,(»)(M+)«(»)n,(»), 


( 2 . 22 ) 


(7=1 

s-hl 


(crm := i:n,(»)*sgti( 


duJr 


M^y^{e)u^ie), eeT^\c,, = o,i. (2.23) 


< 7=1 

Here na.(6') is the spectral projection introduced in Lemma iv). 

Remark 2.8 From Proposition |4.2| , ii) and condition E6 (if Cq 7 ^ 0 ) it follows that 

L^{T^), k,l En. Therefore, (|2.22| ), ( |2.23| ) and ( p.4|) imply that also ((?^)*'^)^; ^ 
k,l En. 

Theorem A Let d,n > 1, a < —d/2 and assume that the conditions El—E5 and SO—S3 
hold. If Co 0, then we assume also that E6 holds. Then 

i) the convergence (\2.^ holds and also holds. 

ii) The limit measure /ioo is a Gaussian translation-invariant measure on Ha ■ 

Hi) The characteristic functional of fioo is the Gaussian 

/ioo(d') = exp{-^Qoo(^, ^)}, ^eV, 

where Qoo is the quadratic form defined in ( \1.13i) . 

iv) The measure /ioo is invariant, i.e. [U{t)]*fioo = Loo , t G IR. 

Remarks 2.9 i) In the case d = n = 1 we have B{9) = 1, and formulas ( 2.14 )- (|2.19|) have 
been obtained in p.l39]. 

ii) The uniform Rosenblatt mixing condition also suffices, together with a higher 
power > 2 in the bound (|1.9|) : there exists h > 0 such that 

e(|mo(x)P+‘^ + |no(a:)P+‘^) < C < oo. 


kl 
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r-\-oo 

Then (|2.13|) requires a modification; / r'^~^a^{r)dr < oo, where p = min(5/(2 + 5), 1/2), 

J 0 

where a{r) is the Rosenblatt mixing coefficient defined as in (|2.11|) but without po{B) in 


the denominator. Under these modifications, the statements of Theorem A and their proofs 
remain essentially unchanged. 


in) The arguments with condition E5 in Lemmas |7.1| , ^^demonstrate that the condition 
could be considerably weakened. Namely, it suffices to assume 

E5’ If for some k ^ I we have either uJkiO) + 0 Ji{ 6 ) = const+ or uikid) — uJi{ 6 ) = const, 
with const± ^ 0, then either pfiiO) - uJk{0)uJi{9)pfi{e) = 0, uJkiO)?^) + ^i{9)Pki{9) = 0 
or p])^{e) + uJk{e)^i{e)PH{ 9 ) = o, uJk{ 9 )pfi{e) - = O. Here 


pSW - (b-(WS(«)B(»))^,, eeT^, k,ten, *,i = o,l, 


(2.24) 


q±{6) are Fourier transforms of covariance matrices q±{z). 

The assertions i)-iii) of Theorem A follow from Propositions |2.10| and |2.11| below. 


Proposition 2.10 The family of the measures {pt, t G IR} is weakly compact in TLa with 
any a < —d/2, and the bounds supE||f/(t)yo||Q < oo hold. 

t>o 


Proposition 2.11 For every T G "P the convergence holds. 


Proposition 2.101 (Proposition |2.11|) provides the existence (resp. the uniqueness) of the limit 
measure poo ■ They are proved in Sections 5 and 7-8, respectively. 

Theorem A iv) follows from (|2.8| ) since the group U{t) is continuous in Tda by Proposi¬ 
tion a). 


2.5 Examples 

Let us give the examples of Eqn ( |1.2|) and measures po which satisfy all our conditions 
E1-E6 and S0-S3, respectively. 

2.5.1 Harmonic crystals 

All conditions E1-E5 hold for ID crystal with n = 1 considered in [^. For any d > 1 and 
n = 1 we consider the simple elastic lattice corresponding to the quadratic form 

d 

{Vu,u) = (^\u{x + Ck) — u{x)\‘^ + m'^\u{x)\‘^), m > 0, (2.25) 

xeZ'" *^=1 

where = {Ski, ■ ■ ■, Skd) ■ Then El holds and = uj‘^{0) with 

uj{9) = 2(1 — cosdi) -|- ... -I- 2(1 — cos9d) + rn^. (2.26) 

Hence, V(x) satisfies E2-E4 with C* = 0. In these examples the set Cq is empty, hence 
the condition E6 is unnecessary. Condition E5 holds trivially since n = 1. Therefore, 
Theorem A holds for (|2.25|) if the initial measure po satisfies conditions S0-S3. In these 
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examples B{6) = 1, and then, q^o = = Mq + isgn {de^oo^M^ , where Mq is dehned 

by (p.l8|) - (p.20|) with ^2(0) = uj{9). Then we obtain the following explicit formnlas for . 
Denote by S{z) = Fq_}^^{ijj~‘^{9)) the fnndamental solntion for the operator —A + m^ on the 

lattice Z'^, i.e. (—A + m?)£{x) = for a: G Z'^, and P{x) = —iFn}^ ^ ^ . Then 

u}[e) 

\ctr + ^ + P * ((q-)°' - (q-)'°);, 

(q+)i° - (q+)°^ +P* ((q“)^^ + (-A + m2)(q-)0°)] , 

(q+)i^ + (-A + m^) ((q+)°° + P * ((q-)°^ - (q-)'°))] , 
where * stands for the convolntion of fnnctions. 


glO = 

^OO 


00 

^zoo 


01 

^oo 


= (-A + = 


1 

2 

1 

2 

1 

2 


2.5.2 Gaussian measures 


We consider n = 1 and constrnct Ganssian initial measnres /xq satisfying S0-S3. We will 
dehne in Tia by the correlation fnnctions (x — y) which are zero for i ^ j , while for 
i = 0,1, 

qt{9) := F^^ebtiz)] e L\T^), qt{9)>0. (2.27) 

Then by the Minlos theorem, there exist Borel Ganssian measnres /i± on Tie , a < —d/2 , 
with the correlation fnnctions g^(x — y ), becanse formally we have 

/ \\y\\l^±{dY) = ^(l + |xn“(trg^°(0)+trgi^(0)) = C{a,d) ftT{f^{e) + q^^{e)) d9 < oo. 

xeZ'' Td 


The measnres y± satisfy SO, S2. Let ns take the fnnctions C± ^ C'(Z) snch that 


a(^) 


1 , for ± s > a, 
0 , for ± s < —a. 


(2.28) 


Let ns introdnee (Y_,Y+) as a nnit random fnnetion in probability space (TdaX-Tda, y-Xy+) ■ 
Then Y± are Ganssian independent vectors in Tda ■ Dehne a “two-temperatnre” Borel 
probability measnre po as a distribntion of the random fnnetion 


Fo(a;) = C-(xd)y-(x) + C+(2:<i)fo+(a;). 


(2.29) 


Then correlation fnnctions of po are 


Qo{x,y) = q"/(x - y)C_{xd)C-{yd) + q+{x - y)C+{xd)C+{yd), hj = 0,1, (2.30) 

where x = (xi,..., Xd ), y = (yi,... ,yd) £ , and q± are the correlation fnnctions of 

the measnres y± . The measnre yo satishes S0-S2. Farther, let ns assnme, in addition to 
that 


<(J(z)=0, | 2 |>r„. 


(2.31) 


Then the mixing condition S3 follows with (p(r) = 0, r > tq . For instance, ( p.27| ) and 
( P-31| ) hold if we set q^-^z) = f(zi)f(z 2 ) ■ ... ■ f(zd ), where f{z) = Nq - |z| for |z| < Nq 
and f(z) = 0 for \z\ > Nq with Nq := [ro/\/d] (the integer part). Then by the direct 
calcnlation we obtain f(9) = (1 — cos Ao6 *)/(l — cos 6*), 9 G T^, and (p.27|) holds. 
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2.5.3 Non-Gaussian measures 

Let us choose some odd bounded nonconstant functions /°, G G(]R). Define /Iq as the 
distribution of the random function {f^{Y^{x)), f^{Y^{x))), where {Y^,Y^) is a random 
function with a Gaussian distribution /xq from the previous example. Then S0-S3 hold for 
Pq with corresponding mixing coefficient </3*(r) = 0 for r > . Measure Pq is not Gaussian 

if the functions are bounded and nonconstant. 

3 Application to Second Law 

We apply Theorem A to the case when are the Gibbs measures corresponding to distinct 
positive temperatures T_ 7 ^ . We deduce that for the limit mean energy current joo = 

we have = —G(T+ — T_) with G > 0. Moreover, under the additional 
condition on V we obtain = 0, k = — 1. This means that the mean energy 

current is directed from high to low temperature in accordance with Second Law. 


3.1 Energy current 

3.1.1 Energy current for finite energy solutions 

We derive formally the expression for the energy current of the hnite energy solutions u{x, t) 
(see (|1.5|) ). For the half-space := {a: G : a:?; > 0} we dehne the energy in the region 
Dfc (cf ([17^)) as 

^ Y, {u{x,t),V{x-y)u{y,t))y 

By formal calculation, using Eqn ( |1.2|) we obtain 

= {Hx,t),V{x-y)u{y,ty - {u{x,t),V{x-y)u{y,t)U . (3.1) 

\a:eO=,i/eOfc ) 

Here := \ = {a: G Z'^ : Xk < 0} . Introduce new variables: a: = a:' -|- mck , 

y = y' where x', y' G Z"^ with = ?/(, = 0 , ek = (4i, • • •, hd ), A: = 1,..., d. Then 

we rewrite ( |3.1|) in the form 

^k{t) = ^ XI { X {u{x+mek,t),V{x+mek-y'-pek)u{y'+pek,ty 

x',y' 

- X {u{x'+mek,t),V{x'+mek-y'-pek)u{y'+pek,ty] =^j’'{x',t). 

m>0,p< —1 x' 

Here j^{x',t) stands for the energy current density in the direction ca, : by dehnition, 

: = ^X{ X (ii{x'+ mek,t),V{x'+ mek-y'-pek)u{y'+ pek,ty 

y' m< —l,p>0 

- X {u{x' + mek,t),V{x +mek-y'-pek)u{y' + pek,t)y^, 

m>0,p <—1 

where x',y' G Z'^ with x'k = y'k = 0 ■ 


12 


3.1.2 Limit mean energy current 


Now let u{x, t) be the random solution to ( [1.21) with the initial measure /tq satisfying S0-S3. 
Then the bounds El and ( |5.3| ) (see below) imply for the mathematical expectation: 

Y. iQr)a(3ix + mek,y'+pek)Vaf3ix'-y'+ {m-p)ek) 

y' m< — l,p>0 

- Y {Q]°)ap{x'+ mek,y'+ pek)Vc,p{x'-y'+ {m-p)ek)). 

m>0,p< — l 

Here we omit the summation on repeating indices a, f3 G n . Therefore, from the convergence 
( p..ll| ) it follows that in the limit f —>■ cx) we get 

Ef{x', t)^j^ = Y -y + {m-p)ek)Va(3{x' - y + {m-p)ek) 

y' m< —l,p>0 

- Y -y'+ {m-p)ek)Vc,p{x'-y'+{m-p)ek)). 

m>0,p< —1 

Denote hy x' — y' =: z', m — p =: s and changing the order of the summation in the series 
we get 

jL = -^YY + sek)V^0{z' + sek)s =-\ Y i.Q^^)c.y{z)zkV^p{z) 

sEl} 

= j^q^^)^^^e)dkV^fs{e)de, fc = i,...,d. (3.2) 

3.2 Gibbs measures 

3.2.1 Definition of the Gibbs measures 

Formally Gibbs measures g± are 

1 + (Vmo,Mo)) 

g±{duQ,dvQ) =—e ^ WduQ{x)dvo{x), 

X 

where (3± = , T± > 0 are the corresponding absolute temperatures. We introduce the 

Gibbs measures g± as the Gaussian measures with the correlation matrices dehned by their 
Fourier transform as 

= = <;"(«) = 9±(«) = o. ( 3 . 3 ) 

Let be the Banach space of the vector-valued functions u{x) G IR” with the hnite 

norm 

ii“ii« = XI (1+ 
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Let us fix arbitrary a < —d/2 . Introduce the Gaussian Borel probability measures g^{du ), 
g]-{dv) in spaces iLQ,(Z‘^) with characteristic functionals (/3± = 1/T± ) 


glii/) = J exp{i{u,^/)} glidu) = exp{- 


(V 

2I3± 


} 


giiij) = J exp{i{v,'ip)}gl{dv) =exp {— 


t/eD = Go(Z'^) 01R" 


By the Minlos theorem, |^, the Borel probability measures g^, g]. exist in the spaces 
because formally we have 

/ Mllglidu) = + / Ui{x)ui{x) gl{du) = + trg^°(0) < cx), 


since a < —d/2 and 


trg^‘'(0) = {27r)-^ j ti q'^iO) dO = T±{27r)-^ j trG-'(0) dO < oo. 

'J'd 'J'd 


The last bound is obvious if Cq = 


0 and it follows from condition E6 if Cq 7^ 0 . Similarly, 


/ 



T±n X! |a;p)" < co, a < —d/2. 

xeZ'^ 


Finally, we dehne the Gibbs measures g±{dY) as the Borel probability measures g^{du) x 
g]-{dv) in {Y G Ha ■ Y = {u,v)}. Let go{dY) be a “two-temperature” Borel probability 
measure in Ha that is constructed in section 2.5.2 with fi±{dY) = g±{dY) and Yq be a 
random function with distribution g^. Denote by gt the distribution of U{t)YQ, f G IR. 
Now we assume, in addition, that Cq = 0, i.e. (cf condition E6) 

det V(0) 7 ^ 0, G (3.4) 


Note that in the case of canonical Gibbs measures condition E5’ is fulhlled (see Remark 
in)). Indeed, by (p.3|) we have 


?;;(») = (B-(») 9 ±(«)B(«))„ = n 4 ,. (3.6) 

Hence. p%{0) = (s*(9)® (9)B(9))^^ = 0 tor fc # i, Vi, j . 


Theorem 3.1 Let conditions E1-E4, ^3.4/) hold and a < —d/2 . Then there exists a Gaus¬ 
sian Borel probability measure goo on Ha such that 


He 


9t 


t ^ oo. 


(3.7) 
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Proof Let us denote by Qt{x,y) the covariance matrix of measure gt, t G IR. Note that 
owing to (p.30|) , the matrix Qo{x,y) is a “linear combination” of q±{x — y). Hence, Qo{x,y) 
satishes conditions S0-S2. Therefore, by (p.3|) we have 


\Qo{x,y)\ < Cl + - y)\, x,y e TL 


± 


|) implies 

|«2:“(z)l = T±|F,-yiC-‘(»)|| ~ (1 + kl)-", VAf € IN. 


Condition 


Hence, Lemma b.l| and Proposition |6.1| (with condition E5’ instead of E5, see Remark p.9| 
iiij) are applicable to the correlation matrix Qt{x,y ), since the proof uses only the bounds 


(3.8) 


of covariance (Rl|) , (|4.2|) . These bounds are now provided by the decay (|3.8|) instead of 
mixing condition S3. Hence, Qt{x,y) —>• Qoo{x,y) , as t —>• cx), and the family of measures 
{gt,t G IR} is weakly compact in Ha, ol < —d/2. Hence, the convergence (|3.7|) holds 
because gt are Gaussian measures. □ 


3.2.2 Limit covariance and energy current for the Gibbs measures 

Now we rewrite the limit covariance qoo{d) and the limit mean energy current j^o dehned 
by (|3.2| ) in the case of the initial measure /xq = 9o with = g± dehned above. At hrst, 
by (p.l8|) - (p.20|) and (|0|) we have 




0 n-^{e) 


(3.9) 


_ -t-T T _ T 

where T := ^ , AT := ^ . Then, due to ( p^ and (|2.16|) we obtain 


2 ’ 2 
r00//3\ _ . — 2 


MMh = Ta;,-^(0)4z, = 

M'Mu = = Arsgn(^(9)pp(9)4,, X./en. 

Therefore, from (|2.15|) we get 

g“(0) = TV-\e), C(0)=T, 

f)l ^ 

&) = -£(«) = -* at (3.10) 


Substituting q/^{9) from ( p.lO| ) in the RHS of (|3^ ), we obtain 


= 

JOQ 


G„i(e)sgn(^(e))i.,,-‘(e)BjV(e)4(B.„(«)s'(«)B;y»))'i»'(3'ii) 


2(27r 


'J'd 


Here as before we omit the summation on repeating indices a, (3,^,6 G n. Since B{9) is 
the unitary matrix, we get 

AT ^ r ,duj^ ^duj^ 

h ^ d9. 


f = 

JcC) 


(27r) 


k=l,....d. ( 3 . 12 ) 


■ySrirpd 
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Remark 3.2 i) From ( p.l2|) it follows that = 


AT 

{27iY 


E 


duj^ 


89. 


{0) 


de<{) if T, >T_ 


■y&nrpi 

ii) In some particular cases we have = 0 for k = 1,... ,d — 1: for example, a) if each 

duj 

u!^{9) is even on every variable 6 *i,..., 9d-i , or b) if sgn [-^^{9)] depends only on variable 

9d. For instance, a) and b) hold for the simple elastic lattice, what follows by ( p.26| ). 

Hi) generally is a discontinuous function by (|3.10|) . Therefore, q)^{x) decays as a 
negative power of |a:|. The exponential decay is impossible in contrast with the case of [|6 


4 Bounds for initial covariance 


Definition 4.1 By V = ® IR"^, p > 1 

f{k) = (/i(fc),..., fn{k)) endowed with norm 


n > 1, we denote the space of sequences 

.» = ( E l/WI”)''''. 


The next Proposition reflects the mixing property in the Fourier transforms q± of initial 
correlation functions . Condition S2 implies that {z) are bounded functions. There¬ 
fore, its Fourier transform generally belongs to the Schwartz space of tempered distributions. 


Proposition 4.2 Let conditions S0-S3 hold. Then 

i) For i,j = 0,1, the following bounds hold 

IQo (^)2/)l — C < oc for all x G Z'^, (4.1) 

\Qo\x^y)\ < C<oo forallyeT'^. (4.2) 

Here the constant C does not depend on x, y G . 

ii) (tieC{T<^), z,j=0,l. 

Proof ad i) Conditions SO, S2 and S3 imply by [|I^, Lemma 17.2.3] (or Lemma i) below): 


\QQ{x,y)\<Ceoip^/‘^{\x-y\), x,yeT^. (4.3) 

Hence, ( p.l3|) implies (|0|) : 

E \Qoi^^y)\ < Ceo Y) (4-4) 

yeZ'" 

ad ii) The bound ( |4.3| ) and condition ( |1.8|) imply the following bound: 

\<li(.z)\<Ce„pl\\z\), zeZ", (4.5) 

Hence, from (|2.13|) it follows that qll{z) G 8 , what implies qf G C{T‘^) . □ 


Corollary 4.3 Proposition i) implies, by the Shur lemma, that for any $, T G the 
following bound holds: 


|(Qo(x,l/),4)(a:)0T(p))|<Ci|«l>||z2||xI/||p. 


(4.6) 
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5 Compactness of measures family 


Proposition 2.10 follows from the bound (|5.1|) by the Prokhorov Theorem [^, Lemma II.3.1] 


using the method of ||22|, Theorem XII.5.2], since the embedding Tta C Tip is compact if 

a> {3. 

Lemma 5.1 Let conditions SO, S2, S3 hold and a < —d/2 . Then the following bounds 
hold 


supE||f/(f)yoila < oo. 


t>0 


( 6 . 1 ) 


Proof Dehnition (|2.1|) implies 


E\\Y{-,t)\\l = ^ (1 + |a:n“(trQ°°(a;,x)+trQ^^^(a:,a:)) < cx). (5.2) 




Since a < —d/2 , it remains to prove that 


sup sup ||(5t(x,|/)|| < C < oo. 




(5.3) 


The representation (|10.3|) gives 

Qpx.j/) = E{Y‘{x,t)(l,Y>{v,t))^ Y. T. ef(x-x')Ql‘(x',y')gi‘(y-y') 

x',y'&T^ fcl=0,l 

= (5.4) 

where 

^l{x\t)-=(gf{x-x'),g]^{x-x')), x'e 7 Z.\ ^ = 0,1. 

Note that the Parseval identity, (|10.5|) and condition E6 imply 

ll«>i(..t)ii?. = f = ( 2 i)-''| (ief(9)p + ie;'(9)p) do < c„ < 

Jf'd J'd 


Then Corollary O gives 

\Qt\x,y)\ = \ {Qo{x’,y’),^l.{x ,t) ^{y',t)) \ < C'||<I'^(-,f)||z 2 ||4>^(-,f)||p < Ci < oo, (5.5) 


where the constant Ci does not depend on x, ?/ G , f G IR. 


□ 


6 “Cutting out” of critical spectrum 

We reduce the proof of the convergences (|1.11|) and (|1.12|) by a suitable spectral analysis. 
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6.1 Equicontinuity of covariance 

Obviously, (| 1 . 11 | ) is equivalent to the next proposition. 

Proposition 6.1 Let conditions E1-E6 and S0-S3 hold. Then VI' G V 

^ t^oo. (6.1) 

Let us show that we can restrict ourselves th G 'P®, where T>^ is a subset of functions 
G P with vanishing spectrum in a neighborhood of a critical set C C T'^. For k = 1, ...,n 
dehne the sets 

Zk:={eeT‘^\C,: Ve,Md)=0}. 


Definition 6.2 i) The critical set C := Cq U C* U U” U ^ U” (see E4/ 
a) := {'^ E V : (9) = 0 in a neighborhood of C} . 


The next lemma plays the central role in our arguments although its proof is similar to the 
proofs of Lemmas and since C ^ T'^. 


Lemma 6.3 Let conditions E1-E4 hold. Then mes C = 0 . 


Next, we introduce a norm || ■ \\v in the space V such that i) is dense in V in 
this norm, while ii) the quadratic forms T), t G IR, are equicontinuous in this norm. 

Then it suffices to prove ( |6.1j ) for T G only. 

Definition 6.4 Vy is the space T> endowed with the norm 

||vl/f,.:= lil + \\V-\9)mi9)\^d9, TgP, (6.2) 

'J'd 

which is finite by condition E6. 

The set T>^ is dense in Vy by Lemma |6.3| and condition E6. 

Lemma 6.5 The quadratic forms Qt(\I/, T), t G IR, are equicontinuous in Vy . 

Proof It suffices to prove the uniform bounds 

sup|Qi(T,T)| ^C-IITII^, TgP. (6.3) 

teR 


Dehnition (p.lO|) implies that Qi(T, T) := E\{Y{xfi), T(x))p . Note that 


{Y{x,t),^{x)) = {Yo{x),(!>{x,t)), 

where $(•, t) := F“^[^^*( 6 ')T( 6 ')]. Therefore, Qt(T, T) = Qo (*^('5 ^) 5‘^(■5 ^)) i so 

sup |Qi(T,T)| < (Fsup ||4'(-,t)||p 
teR teR 


by Corollary |4.3| . Finally, by the Parseval identity and (|10.5|) , we get 




(6.4) 

(6.5) 


( 6 . 6 ) 
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6.2 Equicontinuity of characteristic functionals 


The convergence ([1.12|) also it snffices to prove for T G T’® only. This follows from the next 
lemma. 

Lemma 6.6 The characteristic functionals /it(tl'), f G IR, are equicontinuous in Vy ■ 
Proof This lemma follows immediately from Lemma |63| by the Canchy-Schwartz ineqnality: 

|/ii(Ti) - = I / < I 

< 11 (P, - VI/ 2 ) \MdY) < sjj 

= -T2,Ti -T 2 ) < CIlTi -Tsllv. □ 


7 Convergence of covariance for non-critical spectrum 

We prove Proposition |6]5 for T G T>^ . First we split the initial covariance into the following 
matrices 


Q'^i.x^y) 

:= q+(x- 2 /), 

(7.1) 

Q~{x,y) 

:= q"(x-2/)sgn2/rf, 

(7.2) 

Q'^i.x.y) 

:= Qo{x,y)-Q^{x,y)-Q~{x,y). 

(7.3) 


where q"*" = +<?-)) R = 1^^^+ “*?-)• Since the solntion Y{t) to Canchy problem 

( p..2|) admits the representation ( |10.3|) , we have 


Qt{x,y)= {Gt{x-x')Qo{x\y')gJ{y-y')). 

x' 


Next introduce the matrices 

Qtix^y)= {St{x-x')Q‘^{x',y')gf{y-y')y a;,|/GZ'^, f > 0, (7.4) 

for each a = r} , and split Qt{x,y) into three terms: Qt{x,y) = Qf {x,y)+Qf {x,y) + 

Ql{x,y). Below in Lemmas we will prove the convergence of type (0) to a 

limit for each term Q1{x,y). 

7.1 Convergence of 

Lemma 7.1 lim {Qf {x,y),^{x) ® 4/(2/)) = {qi^{x — y),^{x) ® T( 2 /)) for any 4/ G , 
where the matrix is defined by ^2.22^) . 
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Proof At first, let us apply the Fourier transform to the matrix Qf{x,y) dehned by (|7.4| ). 
Then we have Qt{6,6') := Qf{x,y) = Qti0)Q^i0,F)Q'f (—9'), where Q~^{9,9') : = 

y^-e' 

F^^g Q~^{x,y). From (|7T|) it follows that Q~^{9,6') = 6{6 — 6') {2nY Hence, 

y^-e' 

QUO,O') = (2ir)"i(9 - e')g,{0)ii*(e)g;(e). (7.5) 

Here we use that Qj(—6) = Ql{6) by condition E2. Therefore, 

{«^(a:.9).>l'(a:)®>l'(!/)) = (2^)-“((j+(9.9'). $(9) ® ¥(«'))_ 

= (2ir)-''{ft(9)q+(9)<j,-(9),<P(9)®>i-(9)). (7.6) 

Further, we choose certain smooth branches of the functions B{9) and 0 Jk{d) to apply the 
stationary phase arguments which require a smoothness in 9. We choose a hnite partition 
of unity 


M 


1^9m(9) = l, 9esupp», 


(7.7) 


m=l 


where Qm are nonnegative functions from (7“ (T'^) and vanish in a neighborhood of the set 
C dehned in Dehnition |6.2| , i). Further, using ([7.7|) we rewrite the RHS of ( |7.6|) . Applying 
formulas (|10.6|) , (|10.7|) for Gti9 ), 0^9) > we obtain (see Appendix) 

(«(a:,a),'■'(a:) »'!'(!/)) = (21)-“!: f g.„(9)(B(9)fl,(9)B-(9), $(9) 0 ¥(9)) 99, (7.8) 

^ rpd 

where by Rt{9) we denote the 2n x 2n matrix with the entries (cf ( p.0.13| )): 

M9)ki = lY.{^^^iM9)±u;i{9))t[B*{9){q^{9)TC{9)ct{9)C*{9))B{9) 


± 

sin 


kl 


m{u;Y9)±u;i{9))t [h*( 0)((7(0)q+(0) ± q+(0)(7*(0))H(0)](7.9) 

By Lemma and the compactness arguments, we choose the eigenvalues ujk{9) and the 
matrix B{9) as real-analytic functions inside the supp^^m for every m: we do not mark 
the functions by the index m to not overburden the notations. Now we analyze the Fourier 
integrals with ■ 

At hrst, note that the identitites ujk{9) +uji{9) = const+ or ujk{9) —uji{9) = const, with 
the const± Y 0 ^re impossible by condition E5. Furthermore, the oscillatory integrals with 
ujk{9) ± uji{9) Y const vanish as t —> oo. Hence, only the integrals with ujk{9) — uji{9) = 0 
contribute to the limit, since ujk{9) + uji{9) = 0 would imply uJk{9) = uji{9) = 0 which is 
impossible by E4. We enumerate the eigenvalues 0Jk{9) as in (^A|) . Then if fc, Z G (rg—i, To-], 
we have cos(a;fc — uJi)t = 1 for cr = 1,..., s -|- 1. By formula (|10.15|) with q{9) := q’''(6') and 
o, 0, we get 

{Qt{x,y),'^{x) ® T(j/)) 

_H*(0) + ...,T(0)®T(0))d0 


= / iJ,n(S){B(0)[xu{B-(e)M+(e)B(0)) 

'J'd 

= (27r)-'/(9+(9),>i>(9)0¥(9))99 + ..., 


kU 


(7.10) 


'J'd 
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where Mq{9) is defined in ( |2.1<j|) , ” ..stands for the oscillatory integrals which contain 
cos(a;fc(6*) 9si uji{6))t and sin(ci;fc(6') ± a;;(6*))t with uJkiO) 9si ujiiO) ^ const. The oscillatory 
integrals converge to zero by the Lebesgue-Riemann Theorem since all the integrands in ‘’ 
are snmmable, and V{uiki9) = 0 only on the set of the Lebesgne measnre zero. The 

snmmability follows from Proposition |4.2| , ii) and E6 (if Cq 7 ^ 0) since the matrices B{6) 
are nnitary. The zero measnre follows similarly to ( |2.4|) since ujk{d) ±a;/(6*) ^ const. Lemma 
o is proved. □ 


7.2 Convergence of {x, y) 

Lemma 7.2 lim (QL(a;, ?/), ^(a;) (g) T(?/)) = {q^{x — y), (x) (y)) for any T G , 

t —^OO 

where the matrix is defined in l\2.29i ). 

Proof Step 1 At hrst we apply the Fonrier transform to QT{x,y) dehned by ([7.4|) : 

QT(e,e'):^F,^, QTO,y) = g,(e)Q-{0,e')gT(-0'), (t.h) 

y^-e' 


where Q~{0, 6') := Q~{x, y). Similarly to ( |7.6|) and ( |7.8|) nsing the partition of nnity 

( [f.7|) and also formnlas ( |10.1CI| ) and ([I0.11| ) we obtain 


{QT(X,V), '■'(a:) ® '!'(!/)) = (2jr)-“{(3,-(9. «'). «(«) »»(»')> 
= (27r)-“ '!'(«')). 

m,m' 


(7.12) 


where Rfid, O') is dehned in ( p.0.11| ) with Q{9, O') := Q {0, O'). Second, we have ^^^^(sgn y) = 
i PV( ^ (6*/2) ^’ ^ ^ stands for the Canchy principal part and ?/ G . Hence, 

by o, we obtain 


Q-(0, O') = 5{0 - O') (27r)'^-i i PV( 


tg(^d - ^d)/2 




(7.13) 


Here 0 = (0i,..., Od-i), O' = (0(,..., 0'^_^), 0 = {0, Od), O' = {O', 0'^) G . Note that the 
Fonrier transform of Qf (x, y) is more singnlar than of Qf (x, y) (cf. formnlas (|7.5| ) and 
( ff.lll ), ([7.131) ). Therefore it is of key importance that we can restrict onrselves by the 
fnnctions T G . Fnrther, ( |7.13|) and ( jlU.llj ) with Q{0,0') := Q~{0,0') imply 


Rt{0,0')ki = 5{0-0'){27r)'^-fiFV{ 


:) 


'^g{0d-O'd)/2' 

•^{coso;^;t [M^{0,0'))^^ + smuj^it (^M^{0,0')) (7.14) 


Here 00^0,0') := u;k{0)±u;i{e'), Mt{0,0') ■.= B*{0)^{q{e)TC{0)q-{e)C*{e'))B{0'), 
Mfi{0,0') := B*{6)^(C{0)q~{0) ± q~{0)C*{0')'^B{6'). Let ns analyse the snmmands in the 
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e fci* + e , / + \ — e „ 

— - and sin(a;?,t) = - it sumces 

2 V 2? ^ ^ 

to prove the convergence for arising integrals li (t) resp. Jkiit) with resp. e 

(see Step ^ resp. Step 3). 

Step 2 First, we consider the integrals 41 (t). Let us denote, for simplicity of exposition, 
gm = gm{0), gm’ = gm'{0') and := (^°(6'),4';((6')). Also let us denote by Pki{0,e') 

one of the expressions Brk{d){Mf(9,9'))kiBl^{9') with either + or —, and some i = 1,2, 
r,s en. Then (|7.12|) and (|7.14|) give. 




RHS of (|7.12|) . Since cos = 


J±(f): = {27i)-^^{gr^gm^6{9-9'){27iy-HPV 


ig(9d-9'a)/2 




J'd 


ri 


d9'^y9.{7.15) 
e'={9,e'7i ^ 


The integral with PV in the RHS of ([7.15 ) exists since uji{9') are analytic inside the 
supp5fm'(^0 • Changing variables 9'^ ^ 9^ — 9'^ = ^ in the inner integral in the RHS of 
( [f.l5| ) we obtain 


4l(«) = /g„e‘“>'">‘'S4»)(PV f g„,e 


'J'd 


Ti 


±^M0')t Pki{9,9')^s{9') 

tg(e/2) 


9'={e,ed-0 


d^)d9. (7.16) 


From Dehnition it follows that Ve'jjJi{9') ^ 0 for 9' G supp^f^/ C supp T . Next lemma 
follows from |^, Proposition A.4 i), ii)]. 

Lemma 7.3 Let V q^ioi{9) ^ 0 for 9 G suppi^m' and p{9) G C^{T‘^) . Then 


Pi{9,t) := FV g„,>{9,9,-O^T7r77^Pi9,9d-Od^ 


Ti 


tg(e/2) 


du. 


= 27rz5(^/(6')e=^*"^*(®)V(6')sgn(=F^(6'))+ o(l), t ^+oo, (7.17) 

OOd 


for 9 G supp gm', and 


sup \Pi{9, f)| < oo. 

eeT<i,t&u,i£n 


(7.18) 


Applying Lemma [f.3| to the inner ID integral in (|7.16|) , we obtain as t ^ +cx), 

r\ 

/±(t) = -(27r)-"/(7^(0)(/™i(0)e*‘^-('’^V/(^,^)sgn(T^(0))^(0)^.(0)ci0 + o(l), (7.19) 

rpd 


where a;^/(6', 9) = uJk{9) -LiuJi{9). Let us discuss the limits of the integrals 41 {9) as t —> +CX 0 . 
At hrst, we note that the identities uj'fi{9,9) = const+ or uj'i7i{9,9) = const, with the 
const± 7 ^ 0 are impossible by condition E5. On the other hand, the oscillatory integrals 
with u!^i{9,9) ^ const vanish as t —> cx) owing to Proposition ii), E4, E5, E6 (if 
Co 7^ 0) and the Lebesgue-Riemann theorem (as in Lemma O)- Hence, 


Ikiit) -^0, t^ +00. 


(7.20) 
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Similarly, in the case uj^i{0, 6 *) ^ 0 , we have 0 , as t ^ oo . Therefore, only integrals 

with uj^i{0, 6) = 0 , i.e. k, I G r^] , cr = 1,..., s + 1 (see (^75|) ), contribute to a limit. 
Finally, by (p.l7|) , we get 

= -(27r)“'^y 9m9m'XkiPki{0, 0)sgn (^(6'))t^(6')Ts(6') dO + o(l), t +cx). (7.21) 


Step 3 Now consider the integrals J^(t) of type (|7.15|) with e instead of e*^^*'*. Simi¬ 
larly to (|7.15|) - (|7.20|) , we get 

J+(t): = i27r)-^<^{9m9m'6ie-e')i27rY-HPV ' ^ ~ 

= o(l), t —>■ oo. 


(7.22) 


The same decay as t ^ -|-cx) is valid if we substitute in ( [f.22|) instead of for all 
k,l En except when k,l G (r^_i,rcr] ■ For kj G {r„_i,ro-] we have uJk{d) = 0 Ji{ 6 ) . Hence, 
by the arguments of type ( [f.l5|) -( [fT^ ) and (|7.21|) , we obtain 

Ju(t) = (2i)-^''(9„g™.i(9-90(2^)"-‘iPV—;-2_^e-“S%(9,9').<i'.(«) ®F(9')) 

tg(^d-c^d)/2 

= {2'k)~'^ j 9m9m.'XklPkl{d,d)&g^{^^{d))'^r{d)'^s{d)d9 + o{l), t^+oo. (7.23) 

From (17.201) , (|7.22|) it follows that for any k, I E n as f —> cx), 

{ 9 m 9 m'Si 9 - 9 '){ 2 nY-HFV—-^—^cos{uj^it)pki{ 0 ,e'),'^r{d)^^s{ 0 ')) = o(l), (7.24) 

^S[9d-9d)/ 2 

since the signes in ([7.21|) and ( [f.23|) are oppposite. Similarly, by (|7.20|) and (|7.22|) , we have 

{9m9m'5i9-9'){27iY-^iFV——^—^sm{uj^it)pki{9,9'),^r{9)®^si9')) = o(l), t oo. 

tg[Ud-9^)/2 

(7.25) 

The same relation holds if we substitute uij^i in the LHS of ( |7.25|) instead of ui^i for all 
:ept when fc, / G (ro—i,ro-] 

{2n)-^<^{9m9m'5{9-9'){2nY-\FY- 


k,l En except when k,l E (ro—i,ro-] . At last, using (|7.23|) , we get: 

1 


tg{9d - 9'^/2 


J-— sin(a;fc;t) pki{9, 9'), ^r{9) ® ^s{9')) 


pi. . _ 

{27r)~‘^ {9ra9m'Xki * sgu {-^i(^))pki{9 , 9), '^r{9) ® ^s(6')) "h o(l), t -^CX). (7.26) 


■09^ 

Here (see Steps 1, 2) by Pki{9,9) we denote 

P„{0,O) = B.de){Mi-{0,0))^BU0) = B.de){B-{0)M^{e)B{e))^ 81,(0), 
where Mq {9) is dehned by (|2.19|) . 

Step 4 Now we return to the RHS of ( |7.12|) . Let us substitute ( |7.14|) in (|7.12|) . Then 
by (|7.24|) the summands in the RHS (|7.12|) with coscu^^t tend to zero. Further, by (|7.25|) , 
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( [7.2(j| ) only integrals with sina;^;t , k, I E (ro—i, rg-] , cr = 1,..., s + 1 contribute to a limit. 
Finally, ( fTT^ , (|7-14|) and (|7.22|) - (|7.2q) imply, 


{Qt {x,y),^{x) ® ^{y)) 


'duj^ 




■89^ 

= (27r)“'^ {9m9m' q^i9),^{e) 0 d'(0)) + o(l) 

m,m' 

= {<i^{x-y),'^{x)®'^{y))+ o{l), t^+oo. 


□ 


7.3 Convergence of Ql{x,y) 

Lemma 7.4 lim (QT (x,y),'^(x) '^(y))) = 0 for any ih G . 

t—^oo 

Proof. Step 1 We develop the method [|^ p.l40]. Let us dehne (as in ([6 .41) ) 

$(x',t) := ^ QJ{x — x')'^{x). 
xSZ'^ 

Then using ([7.4|) we have, 

(Q[(x,i/),T(x)0T(|/))= ^ ^ Q^{x\y')^{x\t)^{y\t)= ^ Tt{z'), (7.27) 

x'eZ'^ y'eZ'^ z'eZ'^ 

where 

W)-.= Y. Q^{y' + z\y'W + z\tMy\t). (7.28) 

The estimates (|4.3|), ( |4.5|) and dehnition (|773| ) imply the same estimate for Q''{x,y): 
\Q'^{x,y)\ < — y\). Hence, the Cauchy-Schwartz inequality and (|0|) imply 

WZ)] < Y. WW + Z,-i/)\\\<i(y' + z',t)\\<S>{-i/,t)\ 

B'sZ' 

< |4(j,' + 4()l|-I>(V()l<C.0"(|j'|)||>I'fv-, (7.29) 

y'SlS‘ 


where ||4/||y is defined by (|^) . Hence, ( p.l3| ) and condition E6 imply 

|7^,(^')l < C)®) E 0"(k'|) < C, < oo, (7.30) 

and the series ( 7.27 ) converges uniformly in t. Therefore, it suffices to prove that 

lim tFt{z') = 0 for each z' G 2'^. (7-31) 

t —^OO 

Step 2 Let us prove ( |7.31|) . Condition SI and ([77^) imply that Q^{y' + z\y') = q^{z\y'^ + 
z'd^y'd)^ where 

liin q^{z\y'^ +z^,y'^) =t), for (^, 4) G Z'^. (7.32) 

y^^±oo 
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Hence, Ve > 0 there exists so large that \q^{z'+ z'^,y'^\ < e for \y'^\ > N. 

Respectively, decompose the series (ff.28|) into two series: ^t{,z') = X) S • • • + 

S S • • • • By (|t).6|) and condition E6, the hrst series is estimated by 


Y. Y Y. (7.33) 

s'sZ' 

Note that q^{z', y'^ + z'^, y'^ does not depend on y'. Then we can rewrite the second series 
by the Parseval identity as 


\y'd\<^ 


y'£i 


(2^)_2rf+2 ^ f Fy,^smy' + z',t)]Fy,^omy',t)]de. (7.34) 

L/ I ^ . . 


ly'dK^ 


J^d—l 


It remains to prove that the integral in the RHS of (|7.34|) tends to zero as t —> cxd for fixed 
z' G 2“^ and \y'^\ < N. First, let us note that for the integrand in (|7.34|) the following 
uniform bound holds. 


\Fy,^emy' + z',t)]Fy,^gmy>,t)]\<Gie), f > 0, where G(0) G L'(T''-1). (7.35) 

Indeed, rewrite the function Fy,^§[^{y',t)] in the form 

Fy,^gmy',t)] = (271)-^! e-^^‘^y'<^mt)de,= {27r)-^ I G;{9)^(9) dO,. (7.36) 


Therefore, 


F5.-5[>I>(9' + ='.«)lf’y-»i4>(!/'.()ll < C / lie,•(9)11 I»(9)|d9j <Ci / (|6,-(9)f |>I>(9)|=<i9, 


(pl 'J'l 

< C-,y 11(1 + ||V'-‘(»)l|)|»(9)|^<19i := G(9) 

Ti 


(7.37) 


and ( 7.35 ) follows from condition E6. Therefore, it suffices to prove that the integrand in 
the RHS of ( |7.34|) tends to zero as t —> cx) for a.a. fixed 9 G T'^“^. We use the finite 
partition of unity (|7.7|) (remember that T G ) and split the function Fy,^ 0 [^{y',t)] into 
the sum of the integrals: 


Fy,^g[<!>{y\t)]=Y. / <7m(0)e-*'^^^e±“'=(')*a±(0)T(0)d0,, T G I)°. 


(7.38) 


The eigenvalues ujk{9) and the matrices a^{9) are real-analytic functions inside the supp^f^ 
for every m. From Definition i) and conditions E4, E6 it follows that mes {9d G : 
V0^ijJk{9) = 0} = 0 for a.a. fixed 9 G . Hence, the integrals in (|7.38|) vanish as t —> cx) 
by the Lebesgue-Riemann theorem. □ 

Finally, Lemmas Q.II3 and 0 imply the convergence ( |6.1| ) for T G D® . Then ( |6.1| ) 
follows for any T G D by Lemma ^ (see section 6.1). Proposition |6.1| is proved. □ 
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8 Bernstein’s argument 


8.1 Oscillatory integrals and stationary phase method 

To prove (ICT) we evaluate (y(-,t),T) by (U), where 

4.(3:,«) := -F„-y[er(9)'i'(9)l = (air)-*/ d 0 , x e Z". 

'Y'd 

Similarly to ( [7.38|) or ( [7.81) using the partition of unity ( |7.7| ) we get 

m ±,k&n 


( 8 , 1 ) 


( 8 . 2 ) 


where (Ufc(6') and a^(6*) are real-analytic functions inside the supp^f^ for every m. 
Note that <h(t) := is the solution to the ’’conjugate” equation (cf (pTSp , (|1.4|) ) 


$(t) =A'^t), t e IR; 


A'= 


0 -V 
1 0 


(8.3) 


which is obvious in the Fourier transform. Therefore, the solutions R(t) = 

and <F(t) = (d'°(t), <F^(t)) to the equations ( pTS]) and coincide up to order of the 

components. Hence, ^(x,t) has corresponding dispersive properties. 

We will deduce ( |1.12|) by analyzing the propagation of the solution ^(x,t) to Eqn (|8.3|) , 
in different directions x = vt with v G IR^. For this purpose, we apply the stationary 
phase method to the oscillatory integral (|8.2|) along the rays x = vt, t > 0. Then the 
phase becomes {9v ± ujk{d))t, and its stationary points are the solutions to the equations 
n = =FVa;fc(6'). 

Recall that we can restrict ourselves by T G , hence T(6') =0 in the points 9 E T'^ 
with degenerate Hessian Dk{9) (see E4). Therefore, the stationary phase method leads to 
the following two different types of the asymptotic behavior of <F(nt,f) as t —>• cx): 

I. For the velocity v inside the light cone: v = 9xA7uJk{9) , 9 eT'^\C . Then 

^{vtA) = 0{t-^/^). (8.4) 

II. For the velocity v outside the light cone: v ^ ±Va;fc(6') , 9 E T‘^ \ C, k En . Then 


<F(nt,t) = 0{t ^), Vp > 0. 


(8.5) 


The asymptotics of the types I and II allow ns to incorporate the Bernstein-type approach 
developed in [|l| for case d = 1 and in [^, || for continuous Klein-Gordon and wave equations 
for d > 1. We formalize (|8.4|) , ( p3D as follows. 

Lemma 8.1 For any fixed T G the following bounds hold: 

i) sup mx,t)\ < C t-^/\ (8.6) 

ii) For any p > 0 there exist Cpy'j > 0 s.t. 

|<F(x, t)| < Gp(l-f |x|-F |f|)“^, |(r| > yt. (8.7) 
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Proof Consider $(x, t) along each ray x = vt with arbitrary v G IR'^. Snbstitnting to 
we get 


^vt,t) = Y^ ^ G R>°. (8.8) 

"I ±, ken 

This is a snm of oscillatory integrals with the phase fnnctions = 9v ± ujk{(^) and 

the amplitndes (9) which are real-analytic fnnctions of the 9 inside the snpp gm ■ Since 
ujk{d) is real-analytic, each fnnction 0 ^ has no more than a hnite nnmber of stationary 
points 9 G suppgm, solntions to the eqnation v = ^'VuJk{9). The stationary points are 
non-degenerate for 9 G snpp^fm by (|7.7|), Dehnition and E4 since 

det( ^ ^^ ) = ±Dk{9) 7 ^ 0, 9 E snppfif^. (8.9) 


At last, ^(6^) is smooth since T G T*. 
standard stationary phase method |^, 
|x| < ct with any hnite c. 

Fnrther, denote by h := max^max^gj^ 


Therefore, <h(nt,t) = 0{t '^Z^) according to the 
This implies the bonnds ( p. 6 |) in each cone 


18 


max 

eesupp Qr, 


\Vujj, 


Then for |r;| > n the stationary 


points do not exist on the snpp T . Hence, the integration by parts as in yields <h(nt, t) = 
for any p > 0. On the other hand, the integration by parts in ( p.2|) implies similar 
bonnd <h(x,t) = 0[{t/\x\f'^ for any / > 0. Therefore, ( 0 ) follows with any 7 > n . Now 
the bonnds (| 8 .(j|) follow everywhere. □ 


8.2 ‘Rooms - corridors’ partition 

The remaining constrnctions in the proof of (|1.12|) are similar to 0,0- However, the proofs 
are not identical since here we consider non translation-invariant case. 

Let ns introdnce a ‘room-corridor’ partition of the ball {x G Z'’* : |x| < 'ft} , with 7 
from (0. For t > 0 we choose At and G IN. Asymptotical relations between t , At 
and pt are specihed below. Let ns set = A* -|- pt and 

= jht, V = + At, j G Z, Nt = [{'yt)/ht\. (8.10) 

We call the slabs Rl = {x G Z'’* : |x| < Ntht, < Xd < V} the ‘rooms’, C/ = {x G Z'’* ; 
|x| < Ntht, IP < Xd < the ‘corridors’ and Lt = {x E : \x\ > Ntht} the ’tails’. Here 

X = (xi,..., Xd) , At is the width of a room, and pt of a corridor. Let ns denote by Xt fh® 
indicator of the room Rl , of fh® corridor Cl , and pt that of the tail Lt . Then 

+il{x)]+pt{,x) = 1, xeTL^, ( 8 . 11 ) 


where the snm Yl,t stands for 
tat ion: 


Vt-l 

E • 

j=-Nt 


Hence, we get the following Bernstein’s type represen- 


(Fo,$(-,f)) = V [(E),X^<f>(-,f)) + {Yo,eM;t))] + {Yo,vM;t)). 


( 8 . 12 ) 
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Let us define the random variables rj, , It hj 


Then 


ri = {Yo,xm;t)), cl = {Yo,eM;t)), h = {Yo,vM;t)). 

becomes 

{Y„,<S>(;t)) = J2.{rl + 4) + k- 


(8.13) 

(8.14) 


Lemma 8.2 Let SO—S3 hold and T G "D® . The following bounds hold for t > 1 ; 

E\rl\^ < C(T) A,/L Vj, 

EWt? < C'(vl/) pt/t, Vj, 

E\lt\^ < Vp>0. 


(8.15) 

(8.16) 
(8.17) 


Proof (|8.17| ) follows from ( |8.7|) and Proposition |T^ , i). We discuss ( |8.15| ) only, ( ^.16[ ) is 
done in a similar way. Let us express E\rl\^ in the correlation matrices. Dehnition (|8.13|) 
implies 

E\rl\^ = {Qo{x,y),Xt{x)^{x,t) ® Xt{y)^{y,t))- (8-18) 


According to 


, Eqn (|8.18|) implies that 
E\rl\^ < Ct-'^Y.^t{x)\\Qo{x,y)\\ 


x,y 

= Cr'^^Xtix) ^ \\Qo{x,y)\\<C^t/t, 


(8.19) 


where ||Qo(2^; l/)|| stands for the norm of a matrix (fQ^f{x,y)^. Therefore, (|8.19|) follows by 
Proposition ^]2|, i). 


□ 


Now we prove the convergence ( |1.12| ). As was said, we use a version of the Central Limit 
Theorem developed by Ibragimov and Linnik. If Qoo(d',d') = 0, the convergence (|1.12|) is 
obvious. In fact, then. 


\Eexp{z{Yo, <!>{■, t))} - /ioo(^)| = E\ exp{z(Fo, - 1| < ^1(^0, 

< {e\{Yo, $(•, t))|2)= ((Qo(x, y), <h(x, t) ® $( 1 /, t))) = (Q^(T, T))(8.20) 


where Qt(T, T) ^ Qoo(d', d') = 0 , 
may assume that for a given T G , 


CX) . 


Therefore, ( p..l2|) follows from (|6.1|) . Thus, we 


Qoo(^,^)^0. 


( 8 . 21 ) 


Let us choose 0 < 5 < 1 and 





















Lemma 8.3 The following limit holds true: 

+ (y) ) + + y) ^0) t ^ CO. (8.23) 

Proof. Function (p{r) is nonincreasing, hence by (|2.13|) , 

r r 

rY^\r) = dj ds < d J ds<CTp< oo. (8.24) 

0 0 


Furthermore, (|8.22|) implies that ht = pt +Aj 
Then dS^^) follows by (|;2|) and dpI). 


logt ’ 


oo. Therefore, Nt logt. 

ht 

□ 


By the triangle inequality, 

|F;exp{i(yo,$(-,t))}-hoo(^)| < |^exp{i(Fo,<^(-,^))}-^exp{i^^r^}| 

+ |exp{-i^E|r^f}-exp{-lQoo(d^,d/)}| 

+ l^exp{i^^r^}-exp{-i^^E|r^|2}| 

= h + l 2 + h. (8.25) 


We are going to show that all summands Ji, I 2 , I 3 tend to zero as t ^ 00 
Step (i) Eqn (|8.14|) implies 


Ji = |Eexp{i^^rt^}(exp{i^^4+ih} - Ij 
< V Eldl + E\lt\ < V (E|cif)V2 + 


From (Is:^) , (IST^) , (jST^) and 


we obtain that 


Ii < Cpt ^ + CNt{pt/ty^‘^ ^0, t —>• cx). 
Step (ii) By the triangle inequality, 

+) IfilE/O" - E,-E|Wti + i ieCepT - 


= I 21 + I 22 + T 


23 , 


(8.26) 


(8.27) 


(8.28) 


where Qt is a quadratic form with the matrix kernel (fQ]^{x, y)^ . (|6.1|) implies that J 21 —0 . 
As to I 22 , we hrst have that 


I 22 < E \Er 


irll 


(8.29) 


j<i 


The next lemma is a corollary of [IT, Lemma 17.2.3]. 
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Lemma 8.4 Let A, B be the subsets of with the distance dist {A, i3) > r > 0 , and r] 
be random variables on the probability space (Ha, B(Ha), Lo) ■ Moreover, let ^ be measurable 
with respect to the a-algebra cr(A) , rj with respect to the a-algebra ^(B) . Then 
i) \Eiri - EiEri\ < Cab ^^/^(r) if {E\i\^f/^ < a and {E\7]\^Y/^ < b . 
a) \E^ri — E^Et]\ < Cab (^(r) if |^| < a and \r]\ < b a.e. 

We apply Lemma to deduce that /22 ^ 0 as t —> oo. Note that rf = {YQ,Xt^(-,t)) is 
measurable with respect to the cr-algebra u{Rt) ■ The distance between the different rooms 
Rl is greater or equal to pt according to ( |8.1(J|) . Then (^.291 ) and SI, S3 imply, together 
with Lemma [ 8 ^ i) and (|8.15|) , that 

I 22 < CN^p^/\pt), (8.30) 

which vanishes as t ^ cxd because of ( ^.23 ). Finally, it remains to check that J 23 ^0, 
t —>■ 00 . We have 

"t) = E{Y„, t)f = E{J2,(:rl +4)+ if, 

according to ( p.l4| ). Therefore, by the Cauchy-Schwartz inequality, 

V, < \E{Y.rif-E{Y.ri + Y.A + ‘f\ 

< CN,Y,,E\4? + C,{E{Y,/,)f'\N,J2E\c’X + E\lt\ACcE\k\\ ( 8 . 31 ) 

Then (|8.15|) , (|8.29|) and (p.30|) imply 

+ 2Y^^jErirl\ < CN^A./t + C,Ntp^/\pt) < C 2 < oo. 


Now (|8T6|) , (ICT) , (lOIP and yield 

I23 A CiN^pt/t -\- C2Nt(pt/tY^'^ -\- CA ^ 0 , t —> cx). 

So, all terms I 21 , I 22 , I 23 in (|8.28|) tend to zero. Then (^.28|) implies that 


-Qoo(T,T)|^0, t^CX). 


(8.32) 


(8.33) 


Step (Hi) It remains to verify that 


1 , 


h = \Eexp{iJ2yt} - exp{--E 4 ^kt P}l ^ 0 ^ t 


00 . 


Lemma ^.4| ii) yields: 


Nt-l 

|i?exp{i E/t}- n E explirlll 

-Nt 

Nt-l Nt-l 

< |F^exp{iW^*}exp{z ^ } — F^exp{ir)“^‘}F^exp{i ^ rl}\ 

-Nt+l -Nt+l 

Nt-l Nt-l 

+ \E explirf’^^jE exp{i E w}- n E exp{irl}\ 

-Nt-ei -Nt 

Nt-l Nt-l 

< Cifipt)+ \Eexp{i E n ^exp{ir|}|. 

-Nt+l -Nt+l 
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We then apply Lemma |8.4|, ii) recursively and get, according to Lemma 


Nt-l 


\Eexp{iJ2yt} - n Eexp{iri}\ < CNt^{pt) ^0, t ^ oo. 
-Nt 


It remains to check that 

Nt-l ^ 

I n ^exp{irO-exp{--^^E|r|p}| ^ 0, t ^ oo. 

-Nt ^ 


According to the standard statement of the Central Limit Theorem (see, e.g. 
4.7]), it suffices to verify the Lindeberg condition: Vh > 0, 


0. 


Theorem 


at 


oo. 


Here at = and E^f = E{X^f) , where Xa is the indicator of the event |/| > 

Note that (p.33|) and (|8.21|) imply that at —>■ Qoo('h, 7 ^ 0, t —>■ oo. Hence it remains to 
verify that 




V t\ 


0, t ^ oo, for any e: > 0. 


(8.34) 


We check Eqn (^.34|) in Section 9. This will complete the proof of Proposition |2.11 


□ 


9 The Lindeberg condition 


The proof of (|8.34|) can be reduced to the case when for some A > 0 we have that 

|uo( 2 ^)| + |'i^o( 2 ^)| < A < cxo, x G Z'^. (9.1) 


Then the proof of (|8.34|) is reduced to the convergence 


E.bi 




0, t 


CX), 


(9.2) 


by using Chebyshev’s inequality. The general case can be covered by standard cutoff argu¬ 
ments by taking into account that the bound ( 8.15|) for E\rt\^ depends only on cq and (p. 
The last fact is obvious from ( ^.191) and (|4.4|) . We deduce (p^) from 


Theorem 9.1 Let the conditions of Theorem A hold and assume that l\9.]\) is fulfilled. Then 
for any 4/ G T>^ the following bounds hold: 


E\ri\^ < C{^)X^A.i/ffi t > 1. 


(9.3) 


Proof. Step 1 Given four points X ^ X ^ X ^ X ^ G Z*^, set: 

...,x‘^) = E (Lo(3:^^) ® ... ® Yq{x'^)) . Then, similarly to (|8.18|) , Eqns (|9.1|) and (|8.13|) 

imply 

E\rl\^ = {Xt{x^) • ■ ® • • • 8) Hx^,t)). (9.4) 
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Let us analyze the domain of the in the RHS of (|9.4| ). We partition into three 

parts, W 2 , IL 3 and W 4 : 


= U W„ W^ = {x = a;^ a;^) G (Z*^)^ : |x^ - 

i=2 


max \x^ — x^\}. (9.5) 

p=2,3,4 ' IJ V / 


Furthermore, given x = (a;^, a;^, a;^, a;"^) G IW, divide Z'^ into three parts Sj, j = 1,2,3: 
Z'^ = U 5*2 US's , by two hyperplanes orthogonal to the segment [a:^, a:*] and partitioning it 
into three equal segments, where x^ G Si and a:* G S '3 . Denote by x^ , the two remaining 
points with p,q ^ l,i. Set: Ai = {x E Wi : x^ E Si, x^ E S' 3 } , Bi = {x eWi ■. x^, x'^ ^ Si} 
and Ci = {x E Wi x^, x'^ ^ S' 3 } , i = 2, 3,4. Then IW = U U Cj. Dehne the function 
mo‘^^(a;) , x E (Z'^)*^, in the following way: 


m, 




X 


Wi 


- Qo{x^,x^) ® Qo{x\x‘^), X G Ai, 
Mlf\x), X E BiU Ci- 


(9.6) 


This determines m 


(4) 

0 


(x) 


correctly for all quadruples x . Note that 


{Xt{x^) ■■■Xt {x'^)Qo{x^,xP) ® Qoix\ x''), <F(x\ t) ® ® <l>(x^ t)) 

= iXt ix^)xiix^)Qoix^, ^{x\t) ^ <F(x^, t)) {xi(x^)Xt(x^)Qo(x\ x^), <l>(x*, t) ® <h(x'?, t)). 


Each factor here is bounded by C(T) At/t. Similarly to (|8.15|) , this can be deduced from an 
expression of type ( p.lSj) for the factors. Therefore, the proof of ( |973| ) reduces to the proof 
of the bound 


It ■= \ {Xt{x^) ■ ■■Xt{x^Wo\x^,---,x^),^{x^,t) ® ...® <h(x^,f))| < C'(T)A^A^/t^ 

Step 2 Similarly to (|8.19|) , the estimate ( p. 6 |) implies. 


( > 1 . 

(9.7) 

(9.8) 


We estimate mg'^^ using Lemma |8.4| ii). 

Lemma 9.2 For each i = 2,3,4 and all x E Wi the following hound holds: 

|mo'^^(x^,...,x^)| < C'A^(p(|x^ — x*|/3). (9.9) 

Proof. For x G Al* we apply Lemma ii) to IR^"^ ® IR^" -valued random variables f = 
Yo{x^) ® Yo{x^) and p = Yo{x^) ® Yq{x‘^) . Then (|9TD implies the bound for all x E Ai, 

|ino'^^(x)| < C'A‘^ 9 ?(|x^ — x*|/3). (9.10) 


For X E Bi , we apply Lemma ii) to ^ = Ro(^i) and p = Yq^x^) ® Yq{x^) ® Yq{x’‘) . Then 
SO implies a similar bound for all x E Bi, 

|mS,")(x)| = \M^''\x)-EYo{x^)®E(Yoixp)®Yoix<^)®YoiY))\ < CA\{\x^ - Y\/3), (9.11) 
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and the same for all x G C*. □ 

Step 3 It remains to prove the following bonnds for each i = 2, 3,4: 

- x^\/K) (9.12) 

X 

where Xi is an indicator of the set Ihj. In fact, this snm does not depend on i , hence set 
i = 2 in the snmmand: 


Vi{t) <CY. Xt{x^)Xt{x'^)^{\x^ - x^\/?>)Y,Xt{x^)Y.Xt{x^)X2{x). (9.13) 

x^ x"^ 

Now a key observation is that the inner snm in x'^ is 0{\x^ — x‘^\‘^) as X 2 {x) = 0 for 
Ix*^ — x^l > |x^ — x^l. This implies 

Vi{t) < CY,Xt{x^)Y.Xt{x^)^{\x^ - x^\l?>)\x^ - x^\'^Y.Xt{x^)- (9.14) 

x^ x"^ x^ 


Remember that xi{^) is an indicator of the room Rf = {x G : |x| < Ntht, < Xd < V} , 
where Nt = [xt/ht ]. The inner snm in x^ is bonnded as 


J (p(|x^ — x^|/3)|x^ — x^l'^dx^ < C((i) J ^ip{r/3)dr 

|3:^|<7£ 0 

27t 

< Ci{d) snp r'^Lp^^‘^{r/3) f /3) dr, (9.15) 


rG[0,27£] 


where the ‘ snp ’ and the last integral are bonnded by ( ^.24|) and (|2.13) , respectively. There¬ 
fore, (9.12) follows from (9.14). This completes the proof of Theorem |9.1|. □ 


10 Appendix. Dynamics and covariance in Fourier 
space 

Proof of Proposition |2.4| Applying Fonrier transform to (|1.3| ) we obtain 

Y{t) = A{e)Y{t), t G IR, F(0) = To- 

Here we denote 

A{9) = 


( 10 . 1 ) 


0 1 
-v{e) 0 


e G Th 


( 10 . 2 ) 


Note that Y{-,t) G D'{T'^) for f G IR. On the other hand, V{9) is a smooth fnnction by 
El. Therefore, the solntion Y{9,t) of (|10.1|) exists, is nnique and admits the representation 


Y(9,t) = exp (^A{9)t^Yo{9) which becomes the convolntion 


Y{x,t)= Gt{x - x')Yo{x') 
x'YI3 


(10.3) 
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in the coordinate space, where the Green function Gt{z) admits the Fourier representation 
Qt{z) := Fg:^^[exp (^(0)^)] = (27r)“‘^ y exp {A{6)t)d6. (10.4) 

'j'd 


Hence, by the partial integration, Gt{,z) ~ 
\t\ since .4,(6') is the smooth function of 6 
( p.0.3| ) implies Y(t) G Ha ■ 


z 


-p 

rpd 


as \z\ —>■ oo for any p > 0 and bounded 
Therefore, the convolution representation 

□ 


Covariance in Fourier space Note that Gt{d) has a form 




cos fit 
— sin Qt 


sinfff ^ \ 
cosilt j ’ 


( 10 . 5 ) 


where = fl{9) is the Hermitian matrix dehned by (p?2|) . Let C{9) be dehned by (|2.2CI|) 
and I be the identity matrix. Then 

Qt{9) = cosilt I + sin iltC (9). (10.6) 

Denote by Q{x,y) := e(Yq{x) (8) 4o(p)^ , and Qt{x,y) := E(Y{x,t) ® Y{y,t)'^ ■ Hence, 
applying Fourier transform to Qt{x,y) we get 

Qt{9,9') := E,^e,y^_e,Qt{x,y) = gt{e)Qi9,9')gJ{-9'), 

where Q{9,9') := Ex^ 0 ^y^_giQ{x,y). Note that due to condition E2 il^(—9') = il*{9') = 
il{9') and then gj{—9') = gt{9 '), where 

gi{9) := cos Dt / + sin Dt C'*(6'). (10-7) 

Here C* is Hermitian adjoint matrix to C as in ( p.20|) . Then 

Qt{9,9') = cosil{9)t Q{9,9') cosil{9')t + sinil{9)t C{9)Q{9,9')C*{9') sinil{9')t 

+ cosil{9)t Q{9,9')C*{9') sinil{9')t + sinil{9)t C{9)Q{9,9') cosD(0')h(lO.8) 


Now, for simplicity of calculations, we will assume that the set of the ‘crossing’ points 9^ is 
empty, i.e. uJk{9) ^ ooi{9) , Vfc, / G n, and the functions Uk{9) and B{9) are real-analytic. 
For example, this is the case of the simple elastic lattice ( p.25| ). (Otherwise, we need a 
partition of unity (17.71)). Consider the hrst term in the RHS of (110.81). We rewrite it using 
(PP in the form 

cosf2(6')f Q{9,9') cosil{9')t = B(9)(^cosook{9)t A{9,9')ki cosuji{9')tj{9') 


= B{9)^(^{cos{uJk{9)-ui{9'))t + cos{uJk{9)+uJi{9'))t)A{9,9')ki) j^^^_B*{9'), (10.9) 

where A{9,9') := B*{9)Q{9,9')B{9'). Similarly, we can rewrite the remaining three terms 
in the RHS of ( p.0.8|) . Finally, 

Qti9,9') = gt{9)Q{9,9')g:{9') = B{9)Rt{9,9')B*{E), (10.10) 
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where by Rt{0,6') we denote the 2n x 2n matrix with the entries 


J kl 


Rt{9,9')ki := [B*{9)(Q{9,9')TC{9)Q{9,9')C%ff))B{9') 

^ ± 

+ sin(a;fc(0)±o^;(0'))t [b\9)[C{9)Q{9,9')±Q{9,9')C*{9'))b{9')]^^IQ.II) 

In the translation-invariant case Q{x, y) = q{x — y ), Q{9, 9') = 5{9 — 9')q{9) and we get 

Qt{9,9') = 5{9 - 9')B{9)Rt{9)B\9), (10.12) 

where by Rt{9) we denote the 2n x 2n matrix with the entries 


Rt{9)ki = -E{cos(a;fc(0)±a;K0))t [i?*(0)(g(0)TG(0)g(0)(?*(0))i?(0) 


± 

sin I 


kl 


dn(a;fc(0)±a;z(0))t [B\9){C{9)m ±mC*{9))B{9)]^^. (10.13) 

Let us denote by p{9) := B*{9)q{9)B{9). Then by ( ^73]) and (|2.20|) we obtain 

2 ±'^ V PklT(^kPkl^l Pkl^^kPkl^l J 

+ sinK(9)±^,(9))i ( ‘^oT^ioj)}. (10.14) 

V -‘^kPkl ±Pkl^l -^kPkl TPkl^l J ^ 


We enumerate the eigenvalues u!k{9) as in (p.5|) . Then cos(a;fc(6') — uji{9))t = 1 for fc, Z G 
(r<^_i, r„], cr = 1,..., s -Kl, hence 


RtWk, = l[B-(e){m + c{0)mc’{i)))B{e) 


2L 


kl 


+^cos2uk{9)t [B*{9){q{9) - C{9)mC*{9))B{9) 


2 

1 


+-sm2uk{9)t [B*{9)[C{9)m+Qi9)C*{9))B{9) 


J kl 


kl 


(10.15) 
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